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Abstract 

The perturbed Kepler problem is shown to be a bi-Hamiltonian system in spite of the 
fact that the graph of the Hamilton function is not a hypersurface of translation, which is 
against a necessary condition for the existence of the bi-Hamiltonian structure according 
to the Fernandes theorem. In fact, both the initial and perturbed Kepler systems are 
isochronous systems and, therefore, the Fernandes theorem cannot be applied to them. 


1 Introduction 

Over the last few years Magri’s approach to integrability through bi-Hamiltonian structures 
had become one of the powerful methods of integrability of evolution equations applicable in 
studying both finite and infinite dimensional dynamical systems. 

The global, topological obstructions to the existence of a bi-Hamiltonian structure for a 
general completely integrable Hamiltonian system are discussed in [311 m [S]. Some coun¬ 
terexamples were given to show that an existence of a bi-Hamiltonian structure is not always 
satisfied around a Liouville torus for a given Arnold-Liouville system. For instance, Fernandes 
[8] and Olver [13] announced that the perturbed Kepler problem is a completely integrable sys¬ 
tem without a bi-Hamiltonian formulation with respect to non-degenerate compatible Poisson 
structures in contrast with the initial Kepler problem. 

Below we explicitly present a few non-degenerate bi-Hamiltonian formulations of the per¬ 
turbed Kepler problem using the Bogoyavlenskij construction of a continuum of compatible 
Poisson structures for the isochronous Hamiltonian systems [3]. 

A bi-Hamiltonian vector field is one which allows two Hamiltonian formulations 


A = PdH = P'dK . 


( 1 . 1 ) 


Here P and P' are the two compatible Poisson bivectors with vanishing Schouten brackets 


[P,P] = [P,P'] = [P',P'] =0. 


( 1 . 2 ) 


In generic cases bivectors P and P' could be degenerate and Hamiltonians H and K could be 
functionally dependent. However, under an additional assumption one can construct a complete 
sequence of functionally independent first integrals of X [31 El do]. 

The aim of this note is to present a bi-Hamiltonian formulation of the perturbed Kepler 
vector held X dehned by the Hamilton function 


H = 


pI + pI+pI 

2 


1 

r 2r^ ’ 


r = s/, 


(1.3) 


and canonical Poisson bivector 



(1.4) 
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The r~^ correction can be seen as due to an asymmetric mass distribution of the attracting 
body (e.g., the gravitational attraction of the earth on a nearby orbiting artificial satellite) or to 
non-Newtonian perturbation from the theory of general relativity (e.g., the motion of a particle 
in the Schwarzschild spherically symmetric solution of the Einstein equations), see [17) . 


2 Bi-Hamiltonian formulation of the perturbed Kepler 
problem 


Following to we start with a discussion of the bi-Hamiltonian vector fields in terms of the 
action-angle variables. 

According to the classical Arnold-Liouville theory in the action-angle variables Jk and Wfc 
a given Hamiltonian vector field X = PdH has the simple form 


a: : 


Jk = 0, 


Wfe = 


dH 

dJk ’ 


fc = 1,..., n. 


(2.5) 


Here H = 77(Ji,..., Jn) is a Hamilton function and the Poisson bivector P is the canonical 


one 


n 

^ = E 




— A — 
dJk duJk 


( 2 . 6 ) 


The problem of the existence of the action-angle variables in the neighborhood of an orbit, of a 
level set or globally is discussed in [B] , see also the recent review [T^] and references within. We 
will look for a bi-Hamiltonian formulation only in the domain of definition of the action-angle 
variables. 

The vector field X is called non-degenerate or anisochronous if the Kolmogorov condition 
for the Hessian matrix 


det 


dJ^dJk 


^0 


(2.7) 


is met almost everywhere in the given action-angle coordinates. This condition implies that the 
dense subsets of the invariant n-dimensional tori of X are closures of trajectories. 

In [3] Bogoyavlenskij proposed a complete classification of the invariant Poisson struc¬ 
tures for non-degenerate and degenerate Hamiltonian systems, see Theorem 1 and Theorem 8, 
respectively. 

Let us consider one trivial example of the generic Bogoyavlenskij construction for the 
degenerate or isochronous Hamiltonian system. If in the domain of definition of the action- 
angle variables we have some nonzero derivative 


dH 

dJm 


7 ^ 0 , 


we can make the following canonical transformation 

dH 


Jk — Jk 1 
Jm — H , 


^k — ^k ^ -f a iOf] 

dJk 


k ^ m 


Cjm = a ^UJr, 


( 2 . 8 ) 


This canonical transformation does not add new singularities to the initial action-angle variables 
and reduces the Hamiltonian to the simplest form 

H=Jm- 


It allows us to construct bi-Hamiltonian formulation of the initial vector field X with two 
functionally dependent Hamiltonians 

H = Jm and K = g{Jm) , (2.9) 
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but with the non-degenerate second Poisson bivector 

^ aJt au,u 


dJrr 


-1 


d d 
A 


dJrr 


duJr, 


( 2 . 10 ) 


where Pk{Jk) are arbitrary nonzero functions and g{Jm) is such that g' ^ 0. In this case the 
eigenvalues of the corresponding recursion operator N = P'P~^ are integrals of motion only, 
see examples of this type bi-Hamiltonian formulations of the Kepler problem in mm HE]. 

In fact the Bogoyavlenskij theorems are much more powerful, however, in the following 
discussion, we need only this particular case. 


2.1 Action-angle variables for the pertnrbed Kepler problem 

Let us introduce the spherical coordinates r, 9 and (j) 

r = sjP P , 9 = arctan ^ , (j) = arccos 


x'^ P P 


which are the radius, longitude and azimuth, respectively. In order to describe the corresponding 
momenta we can use the so-called Mathieu generating function 


F = PxT sin (j) cos 9 P pyV sin (p sin 9 -|- PzX cos (p , 


so that 


Pr = 


dP 


P0 = 


dF 


P4> = 


OF 

dp 


In the spherical variables Hamiltonian H (11.31) takes the form 


H = 


1 


p: 


,2 , Pe 


pI 


sin^ 9 


1 e 

r 2r'^ ’ 


such that the Hamilton-Jacobi equation H = h has an additive separable solution 

that allows us to introduce the action-angle variables. Let us begin with the definition of two 
additional commuting integrals of motion 




m = p^ 


which are the total angular momentum and the component of the angular momentum along 
the polar axis. Then for iJ = h < 0 we can calculate the action variables 

J(b — 


P = ^ j>P4‘dp = m, 


Jg = 




pgd9 = 


27r 


^2- 


m 

sirP9 


d9 = i — m, 


( 2 . 11 ) 


27r 


A = — (pprdr = — 


27r 


o, 2 e 

2h-\ - 


dr = 




- VFD 
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using the standard integration method [2]. Then the corresponding angle variables can be 
obtained from the Jacobi equations 


dS ,—— 1 + 2r/i 

bJr = = —rpr\/—2h + arccos ■ 


dJr 


V'l + 2h{P + e) 


ds e ( 1 -^ \ icose 

ujg = = , „„ I arccos — , — | + arcsm 


djg v'l + 2/i(£2 + e) 

dS , , TO cot 6 

= —— = (jjg + (j) + arcsm 


— m? 


dJ,p ^ \/ P — m? 

The Hamiltonian H (|1.3I) in these action-angle variables takes the form 

1 


H = -- 


2 {^Jr -\- {Jg 


2 ■ 


( 2 . 12 ) 


(2.13) 


It is easy to prove that the graph of H (12.131) is not a hypersurface of translation in the action 
variables (12.lip [5] in contrast with the initial Kepler problem at e = 0. 

The perturbed Kepler problem at e ^ 0 and the unperturbed Kepler problem at e = 0 are 
degenerate or isochronous systems 


det 


a2i?(Ji,..., J„) 


dJidJk 


= 0 


with well-defined derivative toi H = h < 0 


dH 


According to the Bogoyavlenskij theorem it allows us to get bi-Hamiltonian formulations of 
these systems in the domain of definition of the action-angle variables (|2.1112.121) . 


2.2 Delaunay type variables 

The action coordinates play an important role in classical dynamics because of their adiabatic 
invariance, i.e. invariance under infinitesimally slow perturbations. In the Kepler problem, 
there are few well-known types of orbits and, therefore, there are few types of the action-angle 
variables associated with different orbits. For instance, the Delaunay elements are valid only 
in the domain in phase space where there are the elliptic orbits [5] . On the other hand the two 
families of the Poincare variables, which are the action-angle coordinates in the phase space of 
the Kepler problem, in the neighborhood of horizontal circular motions when eccentricities and 
inclinations are small [7]. There are also Delaunay-similar elements, Poincare-similar elements 
and some other action-angle variables, which are well-defined in the neighborhood of different 
orbits. 

For the perturbed Kepler problem we can also introduce the Delaunay type variables 

Jl — Jcj) 5 '^2 — J(l> T Jg 5 J^ — Jr T , 

^ (2.14) 

Wl = - Wg , a;2 = We- UJr , UJ 3 = Ur . 

vP + e 

Recall that the Delauney variables have a geometrical meaning directly related to the description 
of the orbits and their variations are much more significant for the astronomers than those of 
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Cartesian or spherical variables mum- For e = 0 variables Jk^oJk (12.1411 coincide with the 
classical Delaunay elements (1, g, h, L, G, H): 

J 3 = L = , 013 = 1 = n(t - r), 


J 2 = G = v'a(l - e2), uj2 = g = uj, 


Ji = H = y^a(l — e^) cos i , wi = h = D , 

where n is the mean motion, a is the semimajor axis of the orbit, e is the eccentricity, i is the 
inclination, w is the argument of the perigee, D is the longitude of the ascending node, r is the 
time when the satellite passes through the perigee. 

In the Delaunay type variables the Hamilton function H (12.131) takes the form 

H = ~ (2.15) 

and, therefore, we can construct the bi-Hamiltonian formulation of the perturbed Kepler model 
with the second bivector P' given by (12.101) . For instance, if 

l3i{Ji) = Ji, h{J 2 ) = J 2 , and K = -^, 

second bivector is equal to 

0 0 Ji 0 0 

0 0 0 J2 0 

0 0 0 0 J3 

0 0000 

J2 0 000 

0 -J3 000 

The corresponding recursion operator has three functionally independent eigenvalues which are 
the first integrals. 

In the initial action-angle variables (w^, wg, Jr,Je,J 4 ,) this bivector has a more compli¬ 
cated form 



. 9 d 


0 

0 

0 

-Jl 

0 

0 


p' = 


0 

0 

0 


0 

0 

0 

0 

0 

-77 


0 

0 

0 

0 

-V 

Je 

J(f) 

Jr — yj 

V 

77 

0 

0 

0 

0 

1 

1 

— Je 

0 

0 

0 

0 

0 

J({} 

0 

0 

0 


where 


V = 


VPT7 


£ — Je + Jcf, ■ 


In much the same way we can obtain other bi-Hamiltonian formulations associated with the two 
families of the Poincare type action-angles variables or with other known types of the action- 
angle variables for the perturbed Kepler problem. Recall, for instance, that action variables at 
the first Poicare type coordinate system are equal to 


Z = Je, r = Jr + V£^ + e-J<t,-Je, A = Jr + Vi'^ + e , 


see [7] and references within. 

Using similar action-angle variables for the relativistic Kepler problem we can also obtain 
the non degenerate bi-Hamiltonian formulation in contradiction with the Fernandes statement 

in[H]. 
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3 Conclusion 


Let us duplicate textually the Fernandes theorem from [5]: 

Theorem: A completely integrable Hamiltonian system is bi-Hamiltonian (satisfying (BH)) 
if and only if the graph of the Hamiltonian function is a hypersurface of translation,relative to 
the affine structure determined by the action variables. 

An additional assumption (BH) is that the corresponding recursion operator N = P'P~^ has 
n functionally independent real eigenvalues Ai,..., An. 

This formulation of the theorem is often mentioned in modern literature, see for instance 
[n El n [El na [m US]. Nevertheless, there is some trivial misprint because the author omits 
one more necessary condition of the non degeneracy (I2J|) of the Hamiltonian function, which 
could be found in the assumption ”ND” on the page 5 in [5] and in the proof of the theorem. 

The author forgets about this condition simultaneously in the formulation of the theorem 
and by considering physical examples of the applicability of this theorem. Thus, Fernandes 
[8] proclaims that perturbed Kepler problem does not have a bi-Hamiltonian formulation in 
contrast with the Kepler problem, see page 13 in [8]: 

’’Also we note that for the unperturbed Kepler problem (e = 0) the graph of the Hamiltonian is a 
surface of translation, and so it has a bi-Hamiltonian formulation (on the other hand, one can 
show that the relativistic Kepler problem also does not have a bi-Hamiltonian formulation).” 
Let us repeat that initial and perturbed Kepler systems are degenerate systems and, therefore, 
we can not use the Fernandes theorem for both these systems simultaneously. 

This work was partially supported by RFBR grant 13-01-00061 and SPbU grant 11.38.664.2013. 
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